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Abstract 

A novel interpretation of MOND is presented. For galactic data, in addition to Newtonian 
acceleration, there is an attractive acceleration peaking at Milgrom's parameter ao- The 
peak lies within experimental error where ao = cHo/2ir and Hq is the present-time value 
of the Hubble constant. This peaking may be understood in terms of quantum mechanical 
mixing between Newtonian gravitation and the condensation mechanism. There are five 
pointers towards galaxies being Fermi-Dirac condensates. 

PACS numbers: 04.50. Kd, 98.62.Dm 

1 Introduction 

The objective of this paper is to introduce a model which gives a physical interpretation of 
MOND. The most relevant astrophysical data will be discussed. They contain five clues which 
suggest that galaxies are quantum mechanical condensates. These clues fit neatly together. 

It is well known that the rotation curves of galaxies deviate from Newtonian mechanics. 
The question is why? Famaey and McGaugh have recently provided a review of all aspects 
of the data, with an exhaustive list of references [TJ. Milgrom's approach should need little 
introduction. He proposed in 1983 a modification of Newtonian mechanics [2], [3], where the 
observed total acceleration a depends on the Newtonian acceleration g ^ as 

a = 9N/Kx), (!) 
and x = a / a o] a o is a constant. For small g^, all forms in use for fi go at large radius r to 

a = VdNdo- (2) 
A star with rotational velocity v in equilibrium with centrifugal force satisfies 




here G is the gravitational constant and M the mass of the galaxy; r cancels and 

v A = a GM. (4) 

It is striking that the t> 4 relation to M agrees well with the empirical Tully-Fisher relation 
between observed velocities at the edges of galaxies and the luminosity [I]. McGaugh shows 
that this relation applies over 6 decades of galactic masses from ~ 10 6 to 10 12 M Q after including 
the mass of gas and dust in each galaxy [5J. It is important that in MOND there is only one 
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free parameter do, fitted to all galaxies once a particular form for fi(x) has been chosen. This is 
the order parameter of the condensate. 

Within the last few years, an independent observation of the same phenomenon has appeared 
in globular clusters. These spherical clusters of stars have dimensions of a few parsec (pc), i.e. a 
factor ~ 10 4 smaller than the Milky Way. Scarpa et al. reported initially on two globular clusters 
situated 16-19 kpc from the Milky Way p]. The equilibrium of such clusters is controlled by 
Jeans' Law, which relates the velocity dispersion of stars to their acceleration. Scarpa et al. 
traced the velocity dispersion of 184 stars at large radius, identified as being members of one 
globular cluster (rather than interlopers), and 146 stars in the second cluster. The velocity 
dispersion is maximal at the centre of each cluster. They were able to trace it to a radius r 
twice that where Newtonian acceleration reaches a®. Velocity dispersions deviate rather abruptly 
from Newtonian acceleration as it decreases through a$. Tidal heating by the Milky Way varies 
as r -3 , and is at least one order of magnitude smaller, making its effect negligible. 

Scarpa et al. have made observations of a further 6 globular clusters. Hernandez and Jimenez 
give the algebra relating velocity dispersions of stars to Newtonian acceleration using Jeans' Law 
[7|. Hernandez, Jimenez and Allen report a detailed study of the velocity dispersion profiles of 
all 8 globular clusters [8J. Like Scarpa et al., they conclude that tidal effects are significant only 
at radii larger by factors 2-10 than the radius where MOND flattens the curves. They also show 
that the velocity dispersion a varies with the mass M of the cluster as M -4 within errors; this 
is the expected analogue of the Tully-Fisher relation arising from Jeans' Law. 

This result is independent of luminosity measurements used in interpreting galactic rotation 
curves. In galaxies, the mass M within a particular radius is not easy to determine, and is 
usually taken as the mass where rotation curves flatten out. Further exploration of globular 
clusters is highly desirable. It has not been suggested till now that Cold Dark Matter has any 
significant connection with globular clusters. 

Several forms for the function \x are used today. Effectively they are mappings and the 
question is how to interpret them. Differences between them are illustrated in Fig. 19 of the 
review of Famaey and McGaugh [lj . A particular form due to Milgrom j9] 



will be used here. It has the smoothest variation with x an d is also symmetric about do- At 
present only 8 globular clusters have been considered, but the algebraic form with which they 
are fitted is closely consistent with this equation. Within individual errors (~ 40%), results 
confirm Eq. (5) for accelerations above do, where galaxies are subject to systematic error in the 
variation of integrated mass as a function of radius. 

The standard cosmological model ACDM assumes that condensation of galaxies is a two- 
stage process. Dark Matter, whose nature is not yet identified, condenses gravitationally and 
provides a host in which galaxies condense as a second step. The view adopted here is that 
quantum mechanics replaces Dark Matter via a one-step process. This should not be surprising. 
The gravitational field is carried by gravitons and quantum mechanics plays a crucial role in 
forming Black Holes. 

The parameter ao of Milgrom's model agrees within experimental error with cHq/2tt, where 
H is the present-time value of the Hubble constant. However, the connection is subtle and 
appears only at the final step. A further clue is that the phenomenology of galactic rotation 
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curves appears on a log- log plot. This will be associated here with the Partition function of 
Statistical Mechanics. Later, it will be associated with the idea that quantised bosons or fermions 
play a role in forming a quantum mechanical condensate that requires that observed features 
appear on a log-log plot. The Partition functions for quantised bosons or fermions are 
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where E s are energies of levels in a box and II denotes the product of the numbers of all possible 
quantum levels in the box. The + sign gives Fermi-Dirac statistics and the — sign Bose-Einstein. 
For a refresher course, see Schrodinger's very clear exposition [TO]; logZ — > J2 s E s /kT in the 
classical limit, k is Bolzmann's constant and T the temperature. 
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Figure 1: (a) The full curve shows the result of equation (10); the dashed line Newtonian 
acceleration, and the dotted one a straight line given by the term */ao#; (b) full curve: the peak 
arising from curvature of the full curve in (a); the dotted curve shows a fit discussed in Section 
4; (c) full curve: the energy derived from the full curve of (b). 



Milgrom plots the log of total acceleration against the log of gravitational acceleration for 
systems of different masses [11]. This presentation is followed in figure 1(a). Axes are x = 
log 10 (77v, where is Newtonian acceleration, and y = log 10 a, where a is total acceleration. 
Newtonian acceleration follows the diagonal dashed line. The full acceleration has a bend 
centred at do- The value of ao is close to 10~ 10 m s~ 2 and it simplifies figures and arithmetic 
to adopt this value. The dotted line of Fig. 1(a) follows the asymptotic form ^a D g of the 
acceleration. In earlier work, attention has been mostly to the asymptotic form of the dotted 
line on Fig. 1(a). Astrophysicists have used this to determine the parameter ao- The dashed and 
dotted curves meet at acceleration a . This suggests some form of cross-over between Newtonian 
acceleration and the acceleration ^a g. 

The scheme introduced here is to focus attention on the extra acceleration which is additional 
to Newtonian acceleration on Fig. 1(a), i.e. the difference between the full curve and the dashed 
line for Newtonian acceleration. A detail is that effects of General Relativity in a galaxy the 
size of the Milky Way are at a level of (2 x 10~ )<?jv; this is well below experimental errors. This 
tiny correction is neglected here and simplifies formulae greatly, [Later, it could be re-introduced 
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as a perturbation.] The additional acceleration is taken as d 2 y/dx 2 of Fig. 1(a) and is shown 
enlarged on Fig. 1(b) as the full curve. 

Section 2 works through the algebra of Figs. 1(a) and (b) and shows that galaxies have 
considerable stability. Subsection 2(a) examines other forms of Milgrom's function (J,(x), with 
the objective of assessing errors in the acceleration fitted to Fig. 1(b). 

Section 3 discusses the total extra energy W obtained by integrating the non-Newtonian 
component over x. The way the algebra is done takes into account explicitly the exponential 
in the Partition function, therefore parametrising the total acceleration exactly. However, it is 
necessary to realise that results appear on a logarithmic horizontal scale. Newtonian gravity is 
normally expressed as a function of r; in equations which appear later, it will be transformed 
to a function of x. This is why Milgrom's function \x has been described as a mapping. It is a 
matter of convention that the zero of Newtonian energy is taken at infinity (but ignoring the 
effect of Dark Energy over the Universe as a whole). Gravitational energy E is then negative for 
finite radii r. The extra energy associated with non-Newtonian gravity has the same negative 
sign as Newtonian energy. It is shown on Fig. 1(c). It turns out that the extra energy at 
infinity is — 0.5GM, where G is the gravitational constant and M is taken as the mass of the 
galaxy integrated from zero radius to the radius where the total acceleration is a . This will 
be interpreted as an energy gap arising from the formation of a Fermi-Dirac condensate. The 
factor of 0.5 arises directly from the factor 2 difference in slope of the dotted and dashed lines 
on Fig. 1(a). 

In order to explain the origin of the energy gap, it is desirable to make contact with techniques 
used in quantum mechanics to understand energy levels in nuclei and mesons. These were 
developed independently by Bogoliubov [12] and Valatin [13]. They introduced the idea of 
working in axes which rotate Fig. 1(a) anti-clockwise by the mean angle of the dashed and 
dotted lines with respect to the x axis; the curve is then symmetrical like the upper curve of 
Fig. 3. The solution for total energy is given by the well known Breit-Rabi equation [33]. There 
are reasons given later why the lower curve of Fig. 3 is unlikely to appear in galactic rotation 
curves. 

The quantum mechanical mixing between gravity and the Hubble acceleration is discussed 
in Section 3.1. A consequence of the asymptotic form ^a gN fitted to galactic rotation curves 
is that it produces a weak tail to the Newtonian potential with a logarithmic dependence on 
radius r. 

Condensates are well known in Particle Physics and are reviewed briefly in Section 4. One of 
them, Chiral Symmetry Breaking, appears to proceed in close analogy with the Fermi function 
describing galaxies. The interpretation given to Fig. 1(c) is that the graviton develops an 
effective mass over the curved section of Fig. 1(a). Its parameters are derived with errors of 
±8% of GM. Exactly how this effective mass arises is beyond the scope of the present paper 
and requires the technical expertise of theorists familiar with renormalisation effects and gauge 
invariance, but a suggestion is made how this may go. Section 5 reviews further work needed in 
interpreting galaxies and Section 6 draws conclusions. 
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2 THE MODEL 



From equation (5), 

(g N /a + a /2a) 2 = 1 + (a /2a) 2 (8) 
gjv + a g N = a 2 . (9) 

The cancellation of the two terms (ao/2a) 2 is typical of a dispersion relation where a subtraction 
may be made at any value of a. This emerges as an important point in Section 4. 

It is inconvenient to use axes to the base 10, so they will be replaced in the algebra which 
follows by conversion to ln e . Where necessary in comparing with Fig. 1, the conversion factor 
2.3026 will be restored. From Eq. (9) 



V = ^ e ^g 2 + a g; (10) 
x = ln e g, (11) 

where the suffix is dropped from gjsj • From equation (10), 

dy/dx = {e x + a /2){e x + ao)- 1 (13) 
d 2 y/dx 2 = {a /2)e x {e x + a )- 2 (14) 
d 3 y/dx 3 = {a /2)e x {a - e x )/{e x + a )" 3 ; (15) 

d 3 y/dx 3 goes to zero at g — ao- [There are higher order terms too.] The curvature d 2 y/dx 2 
has a maximum value at x = ao, where (d 2 y/dx 2 ) = 1/8. The starting equation (5) gives the 
smoothest symmetric curve for d 2 y/dx 2 from the alternative forms for Milgrom's /i(x)- 

In real galaxies, perturbations arise from thermal and pressure effects. Some narrow structures 
are observed in rotation curves, for example from bars in large galaxies. However, these may 
be fitted empirically using Poisson's equation for mass structures appearing in the gravitational 
potential. The MOND component 'rides' such structures smoothly, see Figs. 21 and 29 of 
Famaey and McGaugh pQ. An immediate question is why Planck's constant does not appear in 
results for galactic rotation curves. The reason is that galaxies are noisy enough to hide it. 

Using the known value of H , the maximum curvature is expected at ao = (1.113 ± 0.046) x 
IQ-w m s -2^ j e i g ioClo — _g 953 Results are insensitive to the small difference from 10~ 10 . 
McGaugh summarises a large number of papers comparing the Tully-Fisher relation with models 
of galaxy formation [15] . He concludes that gas rich galaxies give the best determination of the 
baryonic masses of galaxies: ao = (1.3 ±0.3) x 10~ 10 m s~ 2 . A slightly lower value 1.22 ±0.33 x 
10~ 10 ms -2 is found by Gentile, Famaey and de Blok [16]. Fig. 1(b) shows d 2 y/dx 2 as a peak 
dW/dx in the acceleration at ao, marked by the dashed line; here W is the energy of the 'extra' 
contribution. The full curve is well approximated by a Gaussian for the acceleration: 

d 2 y/dx 2 = 0.125 exp-[ 7 (x - a ) 2 ] (16) 

with 7 = 1.175, i.e. the Gaussian drops to half-height at 9.6% of the value of x at the dip. The 
conclusion is that galaxies have considerable stability. 
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2.1 Alternative forms for /i 



Other forms for have been used to fit galactic rotation curves. Their effects are rather 
small. Two common examples are (A) /i(x) — + x)> (B) x/ V^- + X 2 - Both forms have 
the effect of moving the peak of the curvature slightly: from x = — 10 to —9.7 for A and to 
—9.85 for B. The height of the peak for A is scaled by a factor 0.78 and the curve becomes 
correspondingly wider, resulting in a tail reaching 0.012 at x = —8.2 and —11.8. The result for 
Fig. 1(c) is that the top of the Fermi function is 4% of GM lower, and the bottom of the curve 
higher by the same amount, but the central part of the Fermi function is unchanged and still 
centred very close to x — —10. Form B produces the converse effect: a higher, narrower peak in 
Fig. 1(b) and a Fermi function beginning closer to the top and finishing closer to the bottom of 
Fig. 1(c), but with its central section undisturbed. 

A detail is that one might consider the possibility that the smooth curve of Fig. 1(a) could 
be replaced by the dashed and dotted lines. This would give rise to a sharp cusp in Fig. 1(b). 
Such cusps are known in Particle Physics, but arise only at the opening of phase space for a 
new reaction channel [TTJ . However, no such threshold exists in galactic phenomena. 

The Gaussian acceleration of Fig. 1(b) can be integrated to determine the corresponding 
energy using the Error function or the closely related Normal Probability Integral as described 
in items 590 and 585 of Tables of Integrals of Dwight [18J. The result is shown by the full curve 
in Fig. 1(c). 




Figure 2: Axes x, y, x' and y'. 



At this point a rotation of axes is introduced with the objective of making the curve of Fig. 
1(a) resemble Fig. 3. This is the Bogoliubov-Valatin transformation mentioned earlier. It has 
the effect of rotating the full curve of Fig. 1(a) through the mean angle of the dashed and dotted 
lines. This is a standard transformation where there is quantum mechanical mixing between 
any two basis states. The rotation is about the point x = —10, y = — 10 where the two straight 
lines of Fig. 1(a) cross: 

x' = (x + 10) cos/3 - (y + 10) sin/3 (17) 
y' = (a; + 10)sin/3 + (?/ + 10)cos/3. (18) 

Substituting the results of Eqs. (10)-(15) gives an exact expression for the curve in x', y' axes. 
It is also convenient to re-express Newtonian acceleration and energy directly in terms of x, with 
the following simple manipulations; from Fig. 1(a) 

g = e x = GM/r 2 (19) 
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VGM/r (20) 
dr dx , . 

- = - T PI) 



The Newtonian potential is 



£ x = -GM/r = -VGMe x/2 . (22) 



3 The relation to a Fermi function 

The acceleration differs in x and x' axes, but the scalar quantity W is independent of the choice 
of axes; this is used as a check on algebraic manipulations. The simplest form for the Fermi 
function is to take 

wW)« [ 1 + ex p(n^)] ( 23 ) 

where Ep is the energy at the centre of the Fermi function and T is an adjustable unknown 
temperature. This gives the full curve shown in Fig. 1(c). 

It is also of interest to explore how stable this solution is. This has been studied in a number 
of different ways; the one giving the clearest results will be presented here. For this purpose, it 
is useful to express the energy distribution including a parabolic dip near temperature T, hence 
allowing for the possibility that a reduced temperature is appropriate there, corresponding to the 
existence of a condensate. Minor details are needed (a) to introduce a small offset Xq arising from 
the term clo^^On, (b) an exponential form factor to cut off the parabolic dip in the denominator. 
The final form is 



W(x') = -0.5/ 



i a I ( 1 - 7V + x ) 

1.0 + exp 1 



1 + a(x' + x ) 2 exp[— r](x' + x o y 



(24) 



W(x') simply adds to the Newtonian potential. The energy range of the Fermi function in Fig. 
1(c) is 

\AE\=0.5GM. (25) 

This is a large effect locally. As mentioned in the Introduction, it is not accidental that the 
numerical result comes out to be 0.5 GM. This can be traced directly to the asymptotic form 
^aog of MOND, i.e. slope 0.5 on Figure 1(a). 

For small a, i.e. a small parabolic dip, the asymmetry of the fit is very small, as displayed 
by the full curves of Figs. 1(a) and (b). At the other extreme, the dashed curve of Fig. 1(b) 
illustrates the effect of a = 0.6, which increases the \ 2 °f the fit to Fig. 1(b) by a factor 10; 
acceleration is less well determined than the Fermi function because it is a gradient term. The 
effect on the fit to Fig. 1(c) is no larger than the line width. The largest changes observed 
in all tests are ±4% of the height of the Fermi function at both top and bottom. One should 
remember that for real galaxies, asymmetries of this variety can be due to the decrease in the 
relevant mass M inside radius r. 
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Figure 3: Sketch of two crossing atomic lines; full lines show them including mixing, dashed 
lines without; E\ and E 2 label the convention for eigenvalues in the case of no mixing. 



3.1 Quantum mechanical mixing between gravity and the condensation mech- 
anism. 

The Breit-Rabi equation was first derived to describe crossing atomic lines in a magnetic field 
[H]. Where they cross, there is quantum mechanical mixing between them if basis states are not 
orthogonal. Fig. 3 sketches curves and eigenvalues E\ and E2 for atomic levels as a function of 
magnetic field. Away from the crossing point, the eigenvalues vary linearly with magnetic field, 
but repel one another where they cross unless the mixing potential V = 0. The interpretation 
put on Milgrom's formula is that there is distinctive mixing. The maximum repulsion between 
levels is at the crossing point. As pointed out in the Introduction, Fig. 1(a) is a similar diagram 
turned 35.8° clockwise. 

Because galaxies are large ensembles of stars, gas and dust, it is necessary to take expectation 
values of quantum mechanical operators. The place where Quantum Mechanics is needed is to 
describe the mixing between the condensation energy and Newtonian energy and the related 
energy eigenvalues. The eigenvalue equation defines the mixing potential V fitted to data: 



</f*>= ' 1 ;> \ =El I 1 . (26) 





iZiitfi + Vtfa = (27) 
#22*2 + ^1 = E* 2 (28) 
(H ll -E){H 22 -E)-V 2 = 0. (29) 

We work in the galactic centre of mass and in axes rotating with it. The rotational energy is 
ignored; if it is included, it simply describes the rotational equilibrium. H n and H 22 are mean 
values summed over the ensemble. 

Data on galactic rotation curves require that they follow the same behaviour over several 
orders of magnitude in mass from large ellipticals to dwarf galaxies. Then 

H u = Ex = -GM/r = -V~GMe x/2 (30) 
H 22 = E 2 = VGMe{x) (31) 
V = VGMW(x'), (32) 



8 



where e refers to an energy possibly given by the Hubble acceleration. The factor \JGM appears 
in all three equations. The common approximation that the gravitational field of a disc galaxy 
is GM/r is also used, where M is the mass inside radius r. This is accurate to < 1% at the large 
values of r of interest. If the variation with r is known experimentally, detailed information on 
the distribution could easily be adopted for each individual galaxy. 
The two solutions for the Breit-Rabi equation are 



E = ± \ ( ^LZ^f + V 2. (33) 



3.2 A long-range logarithmic tail to the Newtonian potential 

If there were no quantum mechanical mixing between Newtonian acceleration and the conden- 
sation mechanism, there would be no structure at the crossing point and, more serious, no 
explanation for the term ^a gN- 

Asymptotically, the total acceleration, taken from MOND, is 

a = a \fg N /ao- (34) 

Since g N = GM/r 2 , 



a ->• ^GMa /r. (35) 
Taking this as —d<j)/dr, where is a potential induced by the mixing, 



(p ->• -y/GMa ln(r/ri). (36) 

Here r\ is the mean radius for this term. The value of r% is necessarily very close to the much 
larger dip caused by dW/dx'. Because do ~ 1CT 10 , this term is very small. However, it does 
explain the asymptotic straight-line at the right-hand edge of Fig. 1(a). 

The interpretation of this term is simple. Mixing between the Newtonian potential and the 
condensation mechanism allows the wave-length of gravitons trapped in the Newtonian potential 
to expand. This lowers the zero-point energy. An analogy is the covalent bond in chemistry. 
In the hydrogen molecule, each of two electrons is attracted to two protons (which themselves 
repel one another). This increases the wave-length of the electrons and reduces their zero-point 
energy. 

Sobouti noticed this long-range logarithmic component and wrote two closely related papers 
[19], [20]. These included effects of General Relativity. This is interesting but not strictly 
necessary at present. The problem was solved by Sobouti as a power series with additional 
empirical terms proportional to 1/r and 1/r 2 . These terms are replaced by our equations. 

There is evidence that the Milky Way has a halo extending to 0(100 kpc) [21] and [22]. 
However, this could be due to gas and dust shared with the local cluster of galaxies. The halo 
can also be explained by Cold Dark Matter. So this halo is presently ambiguous. 

3.3 Solving the Breit-Rabi equation 

The binding energy W is defined with the same sign as E\. From equation (31), 



2E = E 1 + e(x)±J(E 1 -e(x)) 2 + 4W 2 (37) 
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Vgm / 

ln e 10 V 



■e x/2 + e(x 



± J (e*/ 2 - e(x)) 2 + 4W 2 ^) (3* 



dx ln c 10 ^ */(e a: / 2 -e(x)) 2 + 4W 2 



The factor ln e 10 is included to conform with the fact that Fig. 1(a) has been drawn using axes 
which use logarithms to the base 10. For the upper branch of the solution, the minus sign for 
the term involving the square root is required to reproduce the usual Newtonian potential. 



Parameter 


fit 1 


fit 2 


P 


35.78° 


35.78° 


i 


1.852 


1.864 


7) 


> 4.6 


> 4.6 


X 


0.600 


0.572 


a 


0.0 


0.8 


e 









Table 1: Fitted parameters. 



In practice, the Hubble expansion varies over the radius of the Milky Way only by a factor 
~ 2 x 10~ 4 . As a result, e optimises at zero. Fitted parameters are shown in Table 1 for two 
cases. The best fit is with a ~ 0, i.e. negligible curvature of the dip, and the second gives a 
X 2 worse by a factor 10 for purposes of illustration. The first gives the full curves of Figs. 1(b) 
and (c) and the second gives the dashed curve of Fig. 1(b). The parameter 7' is derived from 
7 of Fig. 1(b), but is larger because of the rotation of axes from x,y to x',y', The value of 77 
is expected to be large, so that the curve for W 2 makes a rapid transition to small values near 
x = —12, where the asymptotic form of ^cl^n generates the very weak logarithmic tail of the 
Newtonian potential. Numerically, 77 > 4.6 and its precise value has little effect on the remainder 
of the fit. The value of a is the least well determined parameter; it is what determines T r . 

A Bose-Einstein condensate does not fit the data. In this case, W should vary as T 3//2 at the 
bottom of the dip [23]. For positive x' near x' — 

W(x') = -B(l + \x'\ 3/2 exp-i x' 2 ) (40) 
dW(x')/dx' = 5(1.5|x'|°- 5 -27 , |x / | 5/2 )exp-7 , x /2 ; (41) 

for negative x', the opposite sign of \x'\ is needed in W(x'). The second term in dW(x')/dx' 
wrecks the x' dependence, which fails to fit the observed peak. Near x' = it has a square 
root variation with x' and then, when the second term of Eq. (40) overtakes it, the curve turns 
downwards. This rules out a Bose-Einstein condensate. 



4 An intriguing connection with Particle Physics 

It is useful to review briefly condensates which are known in Particle Physics, since they offer 
analogies which may carry over to galaxies. The idea of a condensate was introduced by Gell- 
Mann and Levy [23] . They were trying to understand the connection between Weak Interations 
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responsible for /3 decay and the Strong Interactions. Parity conservation is violated in Weak 
Interactions but not in Strong Interactions. Feynman and Gell-Mann provided an explanation of 
the Weak Interactions, using an analogy between Electromagnetism and the Weak interactions 
[25] . Today that idea has been expanded to cover all Weak Interactions and is part of the 
Standard Model of Particle Physics, see Section 10 of the Particle Data Book [26J. The Weak 
Interactions contain a Vector current which is exactly conserved and an Axial current which 
is nearly conserved but not exactly. Goldberger and Trieman derived a relation describing the 
axial current based on perturbation theory [27], but Gell-Mann and Levy were concerned that 
their derivation suffers from renormalisation problems. That same problem remains a central 
issue in gravitation. They introduced a Lagrangian where both it and o fields appear in the 
Lagrangian but with different masses; the a decays to tttt and is responsible for most of nuclear 
binding energy. The Lagrangian can be transformed so that the a term replaces the nucleon 
mass. Along these lines, Nambu and Jona-Lasinio introduced the idea that the nucleon mass 
arises from an energy gap like that in the theory of superconductivity [28J. The idea of nuclei 
as liquid drops is the classic example of a condensate. 

The idea of Chiral Symmetry Breaking was introduced by Gasser and Leutwyler [2H] and is 
today used to account for a range of quantitatively accurate results in meson spectroscopy. The 
first step in understanding the precise mechanism was made by Bicudo and Ribiero [30]. Their 
work finds the Bogoliubov-Valatin transformation to be an essential element. Since then, the 
situation has clarified further via Lattice Gauge calculations. The Strong Interaction is mediated 
by gluons obeying Chiral Symmetry. Below a mass of ~ 1 GeV, the gluon acquires an effective 
mass from its interaction with light quarks, which themselves have masses of ~ 4 and 9 MeV. 
Above 1 GeV, there is a cross-over in which Chiral Symmetry is largely restored and the quark 
model reigns supreme, though with small amounts of mixing with meson-meson and/or qqqq 
basis states. A recent paper of Pennington and Wilson gives details including figures showing 
the cross-over at 1 GeV [3~T] . 

An intriguing conjecture is that the graviton-nucleon interaction follows the same pattern 
as Chiral Symmetry Breaking and develops an effective graviton mass over a narrow range of 
accelerations centred on ao and disappearing a little beyond x = —8, where the curvature of 
Fig. 1(a) becomes negligible. For x < logao, there is a smooth transition to the logarithmic tail 
of the potential. 

Let us now return to Figs. 1(b) and (c). What emerges from the range of fits described by 
Table 1 is that Newtonian acceleration is very small near x = compared with that originating 
from the extra acceleration from dW/dx' . This second term dominates by a huge factor 2.1 x 10 4 
at the centre of the curve, where x = log 10 ao- At x = —8, this factor falls to 70 and is still 
falling fast. The conclusion is that the curved part of Fig. 1(a) is the dominant feature and 
Newtonian gravitation in this region is only a small perturbation. 

Consider the effect of this result near the centre of the Fermi function at E = Ep in Fig. 
1(c). If we retain only the dominant terms in W and dW/dx, the results are 



dE/dx ->■ 



Vgm 

Zn e 10 



dW/dx 



(42) 



E -> 



Vgm 

ln e 10 



W. 



(43) 
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Apart from the factor y/ GM /ln e 10, which is a normalisation factor, dE/dx may be interpreted 
as the modulus of a Breit-Wigner resonance with x-dependent width: 

BW = r(x)/2 / w . (44) 

E-E F - iV{x)/2 K ! 

There is a quantum mechanical pole at E — Ep where dW/dx varies fastest. Astrophysical 
situations would be too noisy to reveal the phase variation. However, it is questionable what 
to take as the absolute zero of energy because of possible renormalisation effects, like those 
responsible for the Lamb shift. The energy W starts at zero because of local gauge invariance, 
and its central value is shifted downward by 0.25 GM. This accounts for the form of Eq. (8) of 
Section 2, where two terms (a /2a) 2 cancel. 

Suppose the graviton acquires an effective mass £(x). The Yukawa potential for gravitational 
coupling between two particles of mass m 1 and m,2 is V(r) = — Gmim 2 e~^ x ^ r / r , where the expo- 
nential accounts for the graviton mass. From the observed variation of dW/dx, the exponential 
falls to half-height for a displacement Ax = ±0.48 on Fig. 1(c). 

In subsection 2.1, the dependence of the fit on alternative forms of Milgrom's \i function was 
tested. Although acceleration curves change significantly, the Fermi function remains rather 
stable. It changes significantly only at the ends of the range x = —8 to —12. In Section 3, 
the stability of the fit using the empirical form of W in Eq. (24) has also been studied. The 
conclusion is again that the shape of the Fermi function in Fig. 1(c) is very stable against these 
changes. 

The conclusion from these results is that the central part of the Fermi function is stable, 
but can be perturbed at the edges. In superconductors, a coherence length was introduced by 
Pippard to account for the effects of defects beyond an experimentally observed range [32]. It 
appears that galaxies behave similarly. 

An attempt was made by t'Hooft and Veltman to construct a renormalisation scheme for 
gravitation, but was inconclusive [33J; that problem persists today. It is obviously of interest 
to re-examine this problem in the extreme infra-red where gravitons interact with nucleons to 
form galaxies. One of the points for investigation is that an effective mass for the graviton 
introduces a separation between the s, t and u channels of the Mandelstam diagram, making 
it resemble Chiral Symmetry Breaking; for massless gravitons, the three channels meet at a 
point. This problem needs to be studied using the Schwinger-Dyson approach along the lines 
of Wilson and Pennington [3T], but is a non-trivial problem. It requires a full understanding 
of how the nucleon acquires its mass from coupling to 3 or more quarks. Presently, there are 
no precise Lattice Gauge calculations of the nucleon mass. There is however an alternative 
approach. Thomas developed a model in 1983 explaining features of deep inelastic scattering 
(i.e. large momentum and energy transfers in electron scattering from nucleons) which provided 
a qualitative explanation of the data in terms of a chiral quark model [31] . This has been followed 
up in a recent paper of Burkardt et al. which improves this work in specific ways [35] . Quoting 
the abstract of this paper: 'The results pave the way for phenomenological applications of pion 
cloud models that are manifestly consistent with the chiral symmetry properties of QCD'. 

Suppose the graviton does acquire an effective mass. A feature of tcN elastic scattering is 
that there is a nucleon pole term at s — M% and an Adler zero at s — u, t — rri^. It would be 
interesting to consider the effect of these perturbations on the behaviour of the nucleon using 
the approach of Burkardt et al. At first sight it appears that the form factor in the gravitational 
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potential will weaken gravity. However, what one must remember is that the wave-lengths of 
gravitons on a galactic scale can be very large, easily larger than a star. The graviton then 
behaves nearly as a plane wave and is a coherent amplitude over a large range. A detail is that 
a massive graviton does develop a scalar component. A coherent amplitude can function like 
Cooper pairs in a superconductor and produce an energy gap. A coherence length like that 
introduced by Pippard can arise in many ways. Supernovae act as major perturbations, heating 
sizable volumes. It is also known that so-called chimneys and wormholes provide channels 
through which currents of dust and gas flow. Furthermore, one of the remarkable features of 
galaxies is that they only grow to a certain size. The largest have masses of order 10 14 M Q . 
Clusters of galaxies do not assemble into one super-galaxy. All of these effects point to a 
coherence length associated with the variation of the Fermi function in Fig. 1(c). 

These remarks are frankly speculative, but deserve further study. There have been papers by 
van Dam [36] and Boulware and Deser [37] which claim no-go theorems preventing the graviton 
developing a mass, but it is not clear whether the present situation evades these theorems. 

Let us return to a simpler issue, the missing lower branch of the Breit-Rabi equation. On 
this branch both W and dW/dx' change sign. This change of sign requires that this branch 
describes an excited state rather than a condensate. Such an excited state is likely to decay on 
a time scale much less than that of galaxies, so it is unlikely that this branch will be observable. 
For those wishing to investigate this branch, the procedure is (a) to fit the upper branch as a 
function of x, (b) rotate to x',y' axes, (c) reverse the sign of y' to reach the lower branch, and 
(d) rotate axes back again to x, y. The best place to search for this branch is near the crossing 
point of Fig. 1(a). 

Presently, Cold Dark Matter (CDM) is used by cosmologists as the source of condensation. 
However, it is not clear whether CDM is cored or cusped. Its radius beyond the peripheries of 
galaxies is not well established. It might also be warm. Recently, Boylan-Kolchin, Bullock and 
Kaplinghat have made a study of the condensation of dwarf galaxies using the standard ACDM 
model [35] • They find that this model performs badly, even when mechanisms such as supernova 
feedback are included. Weinmann et al. reach very similar conclusions [39]. An obvious question 
is what the Dark Matter is that has survived for the age of the Universe without decaying. If it 
did decay significantly, the evolution of galaxies would change as a function of their age. Mond 
and the scheme presented here do not have this problem. 

4.1 The relation to Dark Energy 

If MOND successfully models the formation of galaxies and globular clusters, it raises the 
question of how to interpret Dark Energy. In a de Sitter universe, the Friedmann- Robert son- 
Walker model smoothes out structures using a ACDM function which models the gross features. 
These change over the lifetime of the Universe. In this approach, the Unruh temperature [10] is 
very low, 5 x lO -31 ^, because it is smoothed over the time taken for the Hubble acceleration 
to reach the velocity of light. However, if quantum mechanics governs individual galaxies, there 
will instead be fine structure in Dark Energy. It is logical that steps like Figs. 1(c) do not just 
average out, but instead accumulate over all of this fine structure. The sum total of the fine 
structure is what contributes to the conventional Friedmann-Robertson- Walker model of the 
Universe. This is parametrised via the assumed time dependence of the metric on the Hubble 
acceleration. A similar suggestion along these lines has been advanced by Zhang and Li [H] 
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using ideas based on entropic arguments. The 'present-day' Hubble acceleration is the local 
value and varies over the Universe according to the parametrisation by Dark Energy. 

There is a further argument pointing towards the idea that local fine structure is cumulative. 
Peebles and Nusser argue that galaxies condense more rapidly than the standard ACDM model 
predicts [12] . In particular, they point out that the Local Void contains far fewer galaxies than 
ACDM predicts statistically, while there is an unexpected presence of large galaxies on the 
outskirts of the Local Void. Their Fig. 1 is very persuasive in this respect. Only 3 galaxies 
are observed in the Local Void compared with 19 predicted. The Poisson probability for this 
result is of the order 10~ 5 from ACDM. Peebles and Nusser conclude: 'In short, the general 
sensitivity of galaxies to their environment is not expected in standard ideas. It would help if 
galaxies were more rapidly assembled so that they could then evolve as more nearly isolated 
island universes.' Later Peebles considered an additional empirical term added to the A QCD 
model, but comments that the change requires that Cold Dark Matter is cored rather than the 
expected cusped behaviour [13]. 

A natural explanation is that the Local Void gives no contribution to the Hubble mechanism, 
except for its three galaxies. The total energy E is then higher there. On the periphery of the 
Void, there is gas and dust which can form galaxies. This gas runs down the energy gap to 
enlarge galaxies forming there. It is likely that in the early Universe dwarf galaxies form first; 
globular clusters are known to be very old and may be relics of such small galaxies. Galaxies 
then grow by accumulation of dust and g cts . L ctr ge galaxies are known to spin off small ones in 
collisions. 

5 Further work which is needed 

The model proposed here is precise and open to experimental test. The most important and 
simplest is that if Dark Matter is replaced by this model, it is obviously necessary to redo the 
parametrisation of Dark Energy so that it reproduces smoothly what has been parametrised as 
Dark Matter up to now. This does not necessarily require major modifications. The main point 
is to fit the third peak in the spectrum. This is an exercise requiring cooperation of groups with 
the latest data and techniques at their finger-tips. Despite the clues pointed out here that a 
condensate explains galactic rotation curves, this could fail badly. It would be hardly surprising 
if minor modifications are required in the refit to Dark Energy. 

Angus, Famoey and Buote have considered the possibility that neutrinos may condense within 
the gravitational potentials of galaxies [H] . This is possible, but it is difficult to be quantitative 
at present because only relative masses of neutrinos are known, not absolute values; also there 
is a speculative possibility that sterile neutrinos exist, but that is still an open question. 

The mixing amplitude will affect both Strong and Weak Gravitational Lensing. However, 
unless data are used where a distant quasar transmits light through the periphery of a galaxy, 
where the effect of W(x) is large, the lensing effect will be rather small. The weak logarithmic 
tail of Newtonian gravitation also needs to be taken into account. 

What happens in clusters of galaxies needs detailed, laborious calculations using the formulae 
given here; since there are presently claims that Milgrom's formula does not correctly reproduce 
what happens in clusters, this should be revealing. A multi-body interaction between the Hubble 
mechanism and several galaxies is required, i.e. a generalisation of the Breit-Rabi equation. 
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Quite apart from the attraction to acceleration do, there are also strong tidal effects of the 
variety discussed by Angus et al. [H] and Kroupa [35] . 

The same remarks apply to the Bullet Cluster, which needs to be refitted with the model 
proposed here. The calculation becomes a two-centre problem. The Hubble acceleration (at 
red-shift z ~ 0.2) couples to both galaxies, but they also couple to one another, modifying the 
zero-point energy; individual stars in galaxies communicate with one another not only through 
the Newtonian potential but via their Fermi functions. In the Bullet Cluster, each galaxy 
behaves in this way, but there will be coupling between stars "belonging" to each individual 
galaxy. There may be complex interactions between the two galaxies including resonance effects. 
There are new data on the Bullet Cluster [46] which may offer new interpretations. 

6 Concluding Remarks 

The equations given here allow very little freedom - just the 4% perturbations allowed in the 
Fermi function at the top and bottom of Fig. 1(c). There are five clues which point to galaxies 
and globular clusters being quantum mechanical condensates. Firstly, the phenomena of Fig. 1 
appear on a log-log plot, as expected from Statistical Mechanics; Fig. 1(a) can be interpreted in 
terms of quantum mechanical mixing between two crossing eigenstates. Secondly, Fig. 1(c) is 
fitted naturally by a Fermi function with an energy gap 0.5 GM. Thirdly, the asymptotic form 
of the acceleration in Fig. 1(a) generates a logarithmic tail; this requires a quantum mechanical 
explanation. Fourth, the rotation of axes accomodated by Fig. 3 is just what is expected from 
the Boguliubov transformation for a condensate. A Fermi-Dirac condensate fits the data; a 
Bose-Einstein condensate does not. Fifth, the interpretation of the curve of Fig. 1(b) as an 
energy-dependent Breit-Wigner pole suggests a connection with Chiral Symmetry breaking. It 
suggests that the graviton acquires an effective mass in the vicinity of a . The precise mechanism 
is a challenge to theorists and is presently a conjecture. A final point is that the lower branch of 
the Breit-Rabi equation has the opposite curvature to the upper one; the natural interpretation 
of this result is that it corresponds to an excited state which will decay rapidly and will not 
therefore appear in galactic rotation curves. An unanswered question is why the asymptotic 
form of the acceleration is Ja^g. This is taken from experiment. 




The standard approach to the Universe is the Friedmann-Robertson- Walker model where the 
metric is smoothed over local structures and appears in the metric of General Relativity. My 
suggestion is that galaxies create fine-structure and the Friedmann-Robertson- Walker model 
becomes the sum over these structures. This would explain the agreement between ao and 
cHq/2tt in our part of the Universe. 

Those are the essential points. They do reflect successfully the facts which MOND parametrises 
and they do fit together in a way familiar in Particle Physics in a condensate due to Chiral Sym- 
metry Breaking. 

There are two further small comments. Firstly, there have been many papers on the derivation 
of Newton's law from Entropic arguments. Many assume that one can assign a temperature T 
to the Hubble acceleration and a temperature T% for the gravitational potential by adding them 
as the sum of squares, i.e. the random phase approximation. This gives the result for the total 
acceleration: 





(45) 
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where gn is the Hubble acceleration H^r, and r is the distance to the centre of the gravitating 
object. This is a quite different approach to the one proposed here and leads to a very small 
effect. For the Milky Way, it is a maximum of 2 x 10~ 4 cHo/2ir. 

Secondly, Iorio has calculated that effects of MOND on the perihelia of planets in the Solar 
system are about a factor 10 below present experimental errors [J7]. The logarithmic term 
arising from MOND will be of order a , which is very small and will vary exceedingly slowly 
over the solar system. It has been pointed out by Galianni et al. that it may be feasible to 
detect the MOND effect at acceleration a$ in the solar system near Lagrangian points where the 
accelerations of the Sun, Earth and Moon cancel out [48]. The LISA Pathfinder space mission 
A is planned to visit this area in 2013. If measurements of sufficient accuracy could be made, 
they might confirm or modify MOND as a model of the behaviour of galaxies; secondly, they 
have the important potential to measure the shape of the response function through the region 
where the bend appears in Fig. 1. The alternative approach is the study of globular clusters, 
which provide a different perspective to the rotation curves of galaxies. 
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